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We investigate a dynamical magnetoelectric effect due to a magnetic resonance in helical 
spin structures through the coupling between magnetization and electric polarization via a 
spin current mechanism. We show that the magnon has both the dynamical magnetic moment 
AM" and the electric moment AP" (_L AM"), i.e., a dynamical toroidal moment, under 
external magnetic fields, and thus it is named the toroidalmagnon. The toroidalmagnon exists 
in most conical spin structures owing to the generality of the spin current mechanism. In the 
absorption of electromagnetic waves, the toroidalmagnon excitation process generally induces 
a nonreciprocal directional dichroism as a consequence of an interference of the magnetic and 
electric responses. 

KEYWORDS: multiferroics, toroidalmagnon, nonreciprocal directional dichroism, spin current mechanism, 
frustration 



Electromagnons are magnons that accompany an 
electric dipole moment. Electromagnon excitations are 
thus electro-active, i.e., they are induced by the elec- 
tric component of light, and therefore arc classified 
as dynamical magnetoelectric phenomena. 1-3 Electro- 
magnons are observed as resonances in optical spec- 
troscopy at terahertz frequencies in various mate- 
rials, e.g., i?Mn0 3 , 1 ' 4 ' 5 i?Mn 2 5 , 6 Ba 2 Mg 2 Fei 2 22 , 7 
CuFei_ x Ga x 2 , 8 and Ba 2 CoGe 2 7 . 9 ' 10 Their possible 
applications have been discussed.. 3 ' 11 ' 12 

Electro-active magnetic resonance (MR) processes can 
be understood by spin-dependent electric polarizations 
P({Si}). 13 ' 14 Through the coupling between the electric 
component of light E" and P({Si}), a spin structure 
is modulated by E". When such modulation induces a 
single-magnon excitation in an ordered magnet, the ex- 
citation is called the electromagnon. Thus far, the elec- 
tromagnons induced through the spin current, 2 ' 15 ' 16 ex- 
change striction, 5 ' 17 ' 18 and metal ligand hybridization 
mechanisms 19 have been proposed. 

Although most of these mechanisms have restrictions 
with respect to lattice symmetry, the spin current mech- 
anism can generally induce electric polarization in cy- 
cloidal spin configurations, 20 ' 21 irrespective of local lat- 
tice symmetries. 22 Let us consider a helical magnet with 
a cycloidal spin configuration h _L qo, where qo is 
the propagation vector of the magnetic structure and 
h = Si x Sj gives the hclicity vector. Here, Sj and Sj 
are spins on adjacent sites along the qo direction. The 
spin current mechanism induces the electric polarization 
P ~ (h x qo). In ref. 2, it is discussed that the uniform 
rotation of the helicity vector around the qo direction ac- 
companies the dynamical fluctuation of the electric po- 
larization; thus, such magnon excitation is electro-active. 
Since such magnetic excitation is the Nambu-Goldstonc 
mode, the electromagnon should ubiquitously be identi- 
fied in cycloidal spin systems. 

Another novel feature in the presence of dynamical 



magnetoelectric effects is the nonreciprocal directional 
dichroism (NDD) where absorption intensity depends on 
the direction of an electromagnetic wave propagation 
vector. 23 ' 24 Recently, NDD in MR has been observed in 
Ba 2 CoGe 2 07 9 owing to the electro-active magnetic exci- 
tation process via the metal ligand hybridization mech- 
anisms. 19,25 Since the observation of NDD at a magnon 
resonance is a direct evidence of the presence of dynam- 
ical magnetoelectric effects including an electromagnon 
process, the investigation of NDD is important for under- 
standing the magnetoelectric properties of the system. 
The metal ligand hybridization mechanism, however, de- 
pends on the local lattice structure; thus, NDD due to 
such a mechanism should be observed in a limited class 
of materials. On the other hand, in the generic electro- 
magnons induced through the spin current mechanism 
in cycloidal spin structures, NDD is not observed with 
symmetrical restriction, despite the fact that dynamical 
magnetoelectric couplings exist. 2 

In this Letter, we discuss that NDD should ubiqui- 
tously be observed if one applies a static magnetic field 
to helical magnets with spin current couplings. We clar- 
ify the existence of magnetic excitations that accompany 
both magnetic and electric moments, which contribute 
to dynamical magnetoelectric effects at their resonances. 
We discuss the existence of the dynamical toroidal mo- 
ment in magnon excitation as the origin of NDD. 

Let us first review NDD briefly. 9, 19 In the presence of 
dynamical magneto-electric effects, we have 



= (A c 

= (e°< 



(1) 
(2) 



where x mm and x ee are dynamical magnetic and electric 
susceptibility tensors, respectively, and x me and X em are 
dynamical magnetoelectric susceptibility tensors. As an 
example, let us consider a linearly polarized light with 
E w ||a;, H w ||y, and the propagation vector k||z. In the 
presence of off-diagonal magnetoelectric susceptibilities, 



i 



2 J. Phys. Soc. Jpn. 



Letter 



S. Miyahara and N. Furukawa 



the interferences between E% and Hy in and D w 
given by eqs. (1) and (2) are affected by the relative 
sign of E% and Hy and thus by the direction of k since 
cj/ioEP = k x E". For simplicity, we assume isotropic 
and constant e°° and fx°°. From the Maxwell equation, 
the complex refractive index N is approximated as 



N(s k ) 



X 



=x( W )][Mao 



. v mm 



M] 



-«*[x^H + ^H]/2 



(3) 



as a function of the direction of the propagation vector 
Sk = k z /\k z \. The absorption coefficient is proportional 
to ImN(sk)- Therefore, the absorption intensity depends 
on the light propagation direction Sk and their difference 
is proportional to ImiV(+) — ImJV(— ) = Im[x^(w) + 
X%™(uj)}- In general, when we have off-diagonal magneto- 
electric susceptibilities, 



Im[ X ^»+X e ^M]^0, 



(4) 



the system exhibits NDD for an incident light with 
and Hp for a, (3 = x,y, z (a JL (3). 

Using the Kubo formula, the dynamical susceptibilities 

X^M, X 7 tM, X 7 tM, and X 7 rH for 7, T = X, y,Z 



are described as 
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Here, |0) is a ground state with the eigenenergy E , \n) is 
an ?i-th magnetic excitation state with E n , and hu n Q = 
E n — Eq. V is the unit volume per spin and N is the 
number of spins. Magnetization is defined by 



(9) 



and AM 7 is the 7-component of the magnetization fluc- 
tuations from the ground state. We only consider the 
contributions to the electric polarization induced by the 
spin current mechanism, 20 



x (Si x Sj), 



(10) 



which exists in generic systems in the presence of spin- 
orbit couplings. Here, e^ Hy is the unit vector connecting 
the nearest-neighbor (n.n.) spins and Sj, and A is a 
coupling constant that is proportional to spin-orbit cou- 
plings. AP 7 in eqs. (6)-(8) is the fluctuation of P. 

From eqs. (5) and (6), we see that electromagnon ex- 
citation exists if (n|AP 7 |0) ^ for a onc-magnon ex- 
citation state \n), while a conventional MR occurs for 
(n|AA/ 7 |0) 7^ 0. Furthermore, from eqs. (4), (7), and 
(8), we have NDD at an electromagnon resonance if 
Re(0|AM T |n)(n|AP 7 |0) ^ with r _L 7 being satisfied 




Fig. 1. (Color online) (a) Schematic view of a j/z-cycloidal state 
with P||z induced by the spin current mechanism. The spin ro- 
tation around the j/-axis accompanies AP"||x. (b) yz-cycloidal 
conical state with P||z and M||a:. The tilting of the cones around 
the y-axis creates AP"||a; and AM w [|z. 



for an excited state |n). 

We consider a j/z-cycloidal screw magnet, i.e., qo||y 
and h||x, at B ox = 0, as an example. This struc- 
ture accompanies the electric polarization P|jz. Through 
the Nambu-Goldstone mode, which modulates the spin 
structure in such a way that the cycloidal plane rotates 
along the y-axis, the modulation of the spontaneous elec- 
tric polarization AP"||x occurs, as shown in Fig. 1(a). 
Then, such a mode is activated by E"||x. 2 Since the 
modulation of the spin structure does not induce a uni- 
form magnetization, such an electromagnon mode should 
not be magneto-active, i.e., AM = 0. Thus, the electro- 
magnon resonance should not accompany NDD. 

On the other hand, under the external magnetic field 
B ex ||x, spins give a conical structure, as depicted in 
Fig. 1(b). Here, we have both the electric polariza- 
tion P||z and the magnetization M||x in the ground 
state, and thus have a static toroidal moment defined 
by T = P x M. 26 If we consider an excited state \n) 
where the conical structure is dynamically tilted along 
the y-axis, 15 both AP"||x and AM"p exist coher- 
ently. When the system is excited by an incident light 
with E w ||x and H w ||z, the tilting motion of the con- 
ical structure of the spins around the j/-axis exhibits 
Re(0\AM z \n)(n\AP x \0) ^ and thus NDD. In a semi- 
classical picture, NDD is observed at a resonance ow- 
ing to a magnon with the dynamical toroidal moment 
T d = AP" x AM", i.e., a toroidalmagnon, irrespective 
of lattice dimensions or symmetries. Note that the NDD 
effect is linearly proportional to AP. Thus, in the weak 
coupling region, it may be detected more easily than the 
electromagnon absorption itself, which is proportional to 
(AP) 2 . 

To obtain microscopic details, we investigate a frus- 
trated (quasi-)one-dimensional Heisenberg model, 



H 



Ji ^2 S ; ■ Sj + J2 S i ■ Sj 



(11) 



as a canonical spin model for hclimagnets. Here, is 
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a spin operator on the i-th site. We assume the ferro- 
magnetic n.n. spin exchange interaction J\ and antifer- 
romagnetic next-nearest-neighbor (n.n.n.) interaction J 2 , 
i.e., J\ < and J 2 > 0. In the absence of the external 
magnetic field B ex , the helical spin state is the ground 
state in the parameter range J 2 /|./i| > 0.25 for clas- 
sical spins owing to the frustration. Hereafter, we con- 
sider qo||y as the propagation vector, i.e., the helical spin 
structure along the y-direction has a 27r/g period, where 
cos go = — Ji/4J 2 . When B cx is applied, a conical spin 
state emerges. The periodicity of the helical structure 
does not depend on B cx . 27 We apply the lowest order 
spinwave expansion from the ground state. 

First, we consider a helical magnet with B ox _L qo. For 
B ex ||:r, we have yz-cycloidal and conical spin structures 
with the helicity vector h = x Sj+i||x. The ground 
state spin structure is described as 

Si = (Scosd, S sm8sm(q ri), S sin 6* cos(gorj)), (12) 

where the cone angle 8 is defined as 



Table I. Non-zero coefficients f° (q) for r = x,y, z and a = m,e 
in a yz-cycloidal conical spin structure with qo|[]/ an d B cx ||h||a:. 

AM y f?(±qo) = i[(l±cos8)c qQ + (1 T cos 0)s M ]/4 
AM, / z m (±go) = ±[(1 ± cos 9)c qo + (1 T cos 6>)s 90 ]/4 



2Scos8 



Ji(l - cos g ) 
• it/ 2 for B cx 



(13) 



J 2 (l - cos2g ) ' 

Note that 8 -4 ir/2 for B cx — > such that we recover 
a pure helical structure. We have P||z and M||a; in the 
ground state. 

Following the procedure indicated in ref. 27, we use 
the linear spin wave theory in cq. (11). We define the 
local coordinates where the ^-direction is aligned to the 
equilibrium spin directions: 

Sf = -fif sinfl + fif cos (9, (14) 

Sf = (Sf cos 8 + Sf sin 8) sm{q on ) + Sf cosfeon), (15) 

Sf = (Sf cos 8 + Sf sin 8) cos(q r i ) - S? sm(q Q n). (16) 

Hereafter, we choose qo > for the ground state. Then, 
the Hamiltonian (11) can be written in terms of local £, 
77, and £ coordinates by using S*~, and Sf . For spin 
operators, we consider the Holstein-Primakoff approxi- 
mations: 



S 



+ _ 



25a, 



5; = \/2Sa|, S^ = S — a\c 



(17) 



By the Fourier transformations of aj and a|, the linear 
spinwave Hamiltonian in the q space is written as 



n = E 



where 



SJ2 [2A q a\a q 



B q (a\a 



t t 

1 n 1 



, (18) 



A„ 



— Ji{cos qo + cos 8 sin go sin g 

— — [(1 + cos 2 8) cosg + sin 2 8] cosg} 

— J2{cos2go + cos # sin 2go sin2g 
-i[(l + cos 2 8) cos2g + sin 2 8] cos2g}, 

-j^- sin 2 0(1 — cos go) cos q 

+ ^sin 2 6»(l - cos2g )cos2g. (19) 



AP* 



f£(±qo) = -{c qo -s, )sin 2 6> sing /2 



Here, the magnon wave vector q is defined in an extended 
Brillouin zone to clarify distinctions for the magnon 
modes. The Hamiltonian (18) can be diagonalized by the 
standard Bogoliubov transformation, 

a q = c q a\ + s q a_ qi a q = a q a'_ q + c q a q , (20) 

where the coefficients c q (c q = c_ g ) and s q (s q 
are given by 



2B a 



45 2 - C 2 



4B 2 - C 2 



(21) 



for C q = A q + A^ q - [(A q + A_ ? ) 2 - 45 2 ] 1 / 2 . We obtain 
H = 2^ h UJ q a \ot q + const, (22) 



where the spin wave frequencies are 



S [(A q - A_ q ) + yJ(A q +A- q ) 



A- q f - 45 2 



(23) 



For the helical spin state, i.e., B ex = 0, A q = A- q so 
that the magnons at q and — g are degenerate. On the 
other hand, the external magnetic field lifts the degener- 
acy owing to A q 7^ A_ qi i.e., hoj q 7^ hoj^ q for B ex 7^ 0. 

When the spinwave expansion of P contains linear 
terms of magnon creation and annihilation operators a q 
and a q , the electric component of light creates a one- 
elcctromagnon excited state, through the perturbation 
term %' = — E"AP. In the present model, we apply 
eqs. (17) and (20) to P and M. Both the magnetization 
(9) and the electric polarization (10) can be described as 



AP T 



E 
E 



9f x B V2SN [/ T m (g)aJ + / T m *(g)a 9 ] , (24) 



XSV2SN [r T (qH + JTteK] • (25) 



Non-zero coefficients for the the lowest order spinwave 
expansion /"(g) for a — e,m and r = x,y,z are sum- 
marized in Table I. We only have non-zero contributions 
from g = ±go magnons, which correspond to the elec- 
tromagnons discussed in ref. 2. Note that, within the 
magnetic Brillouin zone, q = ±go locates at the T point. 
At J3 ex = 0, these g = ±go Nambu-Goldstone modes are 
gapless and thus do not show resonance in reality. How- 
ever, under B cx 7^ 0, the q = ±go magnon excitation has 
the gap huj± qo due to Zccman energy Thus, the mag- 
netoelectric resonance of the magnon can be observed 
in a finite frequency range where a conventional MR is 
observed. 

Within a linear spinwave treatment, the susceptibili- 
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ties x™ m (w), X^-M; X™ and X^tM are S iven as 

(26) 

where C mm = 25( 5 ^b)Vo, C 100 = 2S* 3 A 2 /e , and C* mc = 
C em = 25 2 A 5/ i B vW/^- Namely, when \f™(q)\ 2 + 0, 
the resonance due to the magnon excitation at uj = u q 
is induced by , while \f^{q)\ 2 ^ gives an electro- 
magnon excited by . For a non degenerate magnon 
H ^ Re [/ T m *(<z)^(<z)] + (r J_ 7 ) ensures NDD 

for electromagnetic waves with H" and -E 1 ^ with an in- 
tensity oc C mo . 

Since |/J(±g )| 2 ^ as shown in Table I, these 
magnons give electric resonances for xIe( w ) at a; = u) qo 
and W- go • This indeed indicates that MRs can be driven 
by E u \\x as absorptions by electromagnons. We also see 
|/ y m (±g )| 2 = |/r(±9o)| 2 ¥= 0; thus, H- J_ x induces MR 
at the q = ±qo magnon energy. 27 These modes show 
NDD for an incident light with E w ||x, H u \\z, and k||y, 
since x£> ±go ) = xST(^± go ) « /r(±9o)/^(±9o) + 
/|*(±9o)/™(±9o) 7^ 0. In this manner, the electro- 
magnon at q = ±q$ has a dynamical toroidal moment 
T d 7^ and thus is identified as a toroidalmagnon. At 
B cx = 0, however, the ±q magnons become degener- 
ate; thus, the cross-correlated effects at the ui = ui qo 
and oJ~ qo peaks merge to cancel themselves: X™ e ( w g ) °z 

fr(qo)FM + /r(-9o)/|(-go_)_= o. 

Other components of susceptibilities are also obtained 
from Table I. Magnons at q = ±qo contribute to the off- 
diagonal magnetic susceptibility x™( w ) = ~X"™( W ) at 
B CK 7^ 0, which give a conventional Faraday rotation. For 
magnetoelectric susceptibilities, Xj/z ( w ) an d X^y{ w ) arc 
also non-zero. However, they do not cause NDD since 
Im[x^(u;) + x™M]=0. 28 

By reversing the direction of the external magnetic 
fields 5 ox -> -5 ox , xE?(w ±9D ) ~> -xEVtw) is ob- 
tained owing to 9 — > 7r — and huj± qo — > hui Tqo . Thus, 
reversing the light propagation direction is equivalent to 
reversing the external magnetic field direction for the ab- 
sorption process. 

Next, we investigate the case B ox ||q ||?/, where a 
proper conical spin state is stabilized, as shown in Fig. 2. 
Note that the ground state does not have a sponta- 
neous electric polarization owing to ejj||h||y such that 
T = P x M = statically. If we consider an excitation 
state where helical structures are uniformly tilted along 
the z-axis, i.e., h is rotated around the z-axis to create 
Ah||x, the electric polarization AP" ~ (ey X Ah) \\z is in- 
duced, which couples to the external electric field || 3. 
Furthermore, in the case of B ex ^ 0, AM"||a; is also seen 
in this excited state. Thus, the mode is accompanied by 
the dynamical toroidal moment T d ^ even in the ab- 
sence of the static T, such that NDD for an incident 
light with E w ||z and H"||x is expected. The tilting of 
the conical structure along the x-axis is symmetry-wise 
equivalent. This excitation state will induce AP"||x and 
AM w ||z, and NDD for an incident light with E w ||x and 
H w ||z is similarly expected. 

For details, we perform a similar calculation for the 




Fig. 2. (Color online) Schematic view of a rcz-proper conical state 
with M[|j/ but P = 0. The tilting of the cones around the z-axis 
creates AP^Hz and AM" III. 



Table II. Non-zero coefficients /"(g) for r = x,y, z and a = m,e 
in an xz-proper conical spin structure with qo||B ox ||h||2/. 

AM, /i n (±<?o) = -i[(l ± cos + (1 =F cos0)s, o ]/4 

AM Z /f (±gp) = ±[(1 ± cos e)c q + (1 T cos e)s qo ]/4 

AP X fx(±Qo) = -(cq -s qg )sm 2 e sin<j /2 

AP Z /|(±(}o) = =F*( c 9o -s qo )sin 2 6 singo/2 



spinwave expansions in a proper conical spin state with 
q ||B cx ||h||2/. We obtain AP and AM within the linear 
combinations of spinwave operators, as listed in Table 
II. We see that magnons with q = ±q are electro-active, 
owing to \f°(±q )\ 2 = |/ 2 c (±9o)| 2 ± 0. Namely, these 
modes are electromagnons that couple to E w _L y, irre- 
spective of B ex . Moreover, they are also magneto-active, 
since |/™(±<7o)| 2 - l/ z m (±<7o)| 2 ? at B cx ^ 0. 

Let us study NDD in this configuration in detail. For 
W\\z and H"||ar, we have Im[x£f(u; ±(R ,) + X^(w±„)] ^ 
0, owing to fr(±qo)m±qo) = P* (±q )f™(±qn) + 0. 
Similarly, we have Im[xf°(u>± qo ) + XxT( w ±?o)] ^ for 
E"||x and H"||z. At B cx = 0, however, the cross- 
correlated effects caused by igo-magnons are canceled 
owing to their degeneracy. Namely, for an incident light 
with k||y, we see NDD for any polarizations of E 1 ^ and 
H", in the presence of B ex ^ 0. Note that this is the 
Faraday configuration k||M such that circularly polar- 
ized modes are the eigenmodes of the Maxwell equation 
in the presence of rotational symmetry along k|jM; thus, 
the so-called magnetochiral dichroism is observed. De- 
tails will be reported elsewhere. As in the cycloidal cone 
case, reversing the light propagation direction is equiv- 
alent to reversing the external magnetic field direction 
for the absorption process: B cx — > —B cx (9 — > it — 9 and 
huj± qo -> fuj Tqa ) induces X™ c ( w ± 9o ) ~> -x'"x°( w t«o) and 

X^>± 9 o)^-X™K<7o)- 

Microscopically, NDD arises owing to the off-diagonal 
matrix elements of both the magnetic and electric dipolc 
operators Re(0| AM 7 |n) (n\ AP T \0) (7 ^ t) [sec eqs. (4), 
(7), and (8)]. In general, the magnetic dipole operator 
is even (odd) under space-inversion / (time-reversal R) 
symmetry operation, whereas the electric dipole operator 
is odd (even) under / (R). Thus, in magnetic structures 
with J and R symmetry invariances, the real parts of 
the off-diagonal matrix elements are zero, i.e, NDD does 
not arise. In fact, the matrix elements in pure helical 
structures vanish owing to the R symmetry invariance, 
and NDD is induced only in conical structures where I 
and R symmetries are broken. 

To summarize, we investigated NDD that arises in 
conical spin structures via the spin-current mccha- 
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nism. NDD is generally induced by the excitation of a 
magnon accompanied by a dynamical toroidal moment, 
i.e, toroidalmagnon. In a conical structure, the fluctua- 
tion of the dynamical toroidal moment is observed with 
cone oscillations, as shown in Figs. 1 and 2. Such modes 
are adiabatically connected to the q = ±q Nambu- 
Goldstone modes, and thus we conclude that a toroidal- 
magnon should ubiquitously exist in conical spin struc- 
tures. Let us finally note that the observation of NDD 
via the magnon excitation in a perovskite helimagnet 
has been reported in ref. 29 quite recently, which is qual- 
itatively explained by our theory. 
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